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Abstract 

In this paper we study the convergence of monotone PI finite element methods for fully 
nonlinear Hamilton Jacobi-Bellman equations with degenerate, isotropic diffusions. The main 
result is strong convergence of the numerical solutions in a weighted Sobolev space L 2 (f/^(fl)) 
to the viscosity solution without assuming uniform parabolicity of the HJB operator, finite 
element methods, degenerate partial differential equations, Hamilton-Jacobi- Bellman equations, 
viscosity solutions. 


1 Introduction 

Hamilton-Jacobi-Bellman (HJB) equations characterise the value functions of optimal control prob¬ 
lems. For a wide range of control problems one can compute optimal control policies from the partial 
derivatives of the value function. 

An important tool in the analysis of HJB equations and their numerical approximations is the 
concept of viscosity solutions. Its definition is based on sign information on function values of 
candidate solutions, leading typically to proofs of L°° convergence of numerical methods, cf. [2], It 
is more difficult to prove convergence in other norms if solely viscosity solutions are used. 

The use of weak solution, familiar from semilinear differential equations, in the context of 
Hamilton-Jacobi-Bellman equations is delicate because often uniqueness cannot be ensured. How¬ 
ever, we believe that combining the notions of viscosity and weak solution is attractive for numerical 
analysis: the former to deal with uniqueness and the later to study convergence of partial deriva¬ 
tives. 

In [10] the uniform convergence of PI finite element approximations to the viscosity solutions 
of isotropic, degenerate parabolic HJB equations was shown. In addition L 2 (H l ) convergence was 
demonstrated, under the assumption that the HJB equation is uniformly parabolic. In this paper 
we remove the assumption of uniform parabolicity and verify that strong convergence in weighted 
L 2 (H^) spaces can be maintained, see Theorem 1 below. Also the condition that the rJ-dimensional 
Lebesgue measure of the boundary of the zero level set of the value function has to vanish is not 
needed anymore. 
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Our approach uses coercivity properties of the HJB operator. An alternative technique to 
control derivative terms is proposed in [11], where HJB equations satisfying Cordes conditions are 
discretised. A general review of the recent advances in the discretesation of fully nonlinear equations 
is given in [7]. 

In Section 2 we introduce the Bellman equation, while in Section 3 the numerical method is 
defined. Section 4 is concerned with uniform convergence to the value function. In Section 5 a 
nonlinear projection operator is constructed and analysed, which preserves positivity and boundary 
conditions of the viscosity solution. Section 6 is concerned with the coercivity of the continuous 
and discrete Bellman operators. The main result of strong convergence in a weighted Sobolev space 
is proved in Section 7. Finally, in Section 8 assumptions of the prior analysis are translated into 
concrete parameter values for the method of artificial viscosity. 

2 Problem statement 

Let fl be a bounded Lipschitz domain in d> 2. Let A be a compact metric space and let 

A -a C{U) x C(H, R d ) x C(U) x C(U), a ^(a a ,b a ,c a , /“) 

be continuous, such that the families of functions {a a } a ^Ai {b a } a eA, {c a } a eA and {f a } a eA are 
equi-continuous. Consider the bounded linear operators 

L a : H 2 {Q) -a L 2 (H), w —a° Aw + b a ■ Vw + c a w , a e A. 

We assume that a° > 0, i.e. that all L° are of degenerate elliptic. Furthermore, suppose that 
pointwise /“ > 0. Then 


sup || (a a ,b a ,c a ,f ) Hc(n)xC(o,R d )xC(f 2 )xC(a) < 00 > (-*-) 

and also sup agj 4 ||-L a ||_ff 2 (f 2 )->T 2 (n) < oo. Let the final-time data vt € C(Q) be non-negative, that is 
vt > 0 on 17, and let vt satisfy homogeneous boundary conditions on dCl. For smooth w, let 

Hw := sup {L a w - /"), 

a 

where the supremum is applied pointwise. The HJB equation considered is 

—dtv + Hv = 0 in Qt ■= (0, T) x H, (2a) 

v = 0 on (0,T) x dQ, (2b) 

v = vt on {T} x n. (2c) 

Definition 1 ([1, 8]) An upper semi-continuous (lower semi-continuous) function v : [0, T] xll —> M. 
is a viscosity subsolution (supersolution) of 

—dtv + Hv = 0 on Hr, (3) 

if for any w € C' 0O (lR. X ]R d ) such that v — w has a strict local maximum (minimum) at (t,x) G 
(0,T) x H with v(t,x) = w(t,x), gives —dtu>(t,x) + Hw(t,x ) < 0, (greater than or equal to 0). 
If v € C([0,T] x H) is both a viscosity subsolution and a supersolution of equation (3), then v is 
called a viscosity solution. 


2 



The viscosity solution of (2) is understood to be a viscosity solution of the PDE (2a), in the 
sense of Definition 1, that satisfies pointwise the boundary conditions (2b) and (2c). Owing to the 
definition of viscosity solutions, v is a continuous function. 


3 The numerical scheme 


Let Vi, i € IN, be a sequence of piecewise linear, conforming, shape-regular finite element spaces 
with nodes yf. Here £ is the index ranging over the nodes of the finite element mesh 7). Let 
V) 0 C Vi be the subspace of functions which satisfy homogeneous Dirichlet conditions on dfi. It is 
convenient to assume that yf € ft for £ < Ni := dim If 0 ; i.e. the index £ first ranges over internal 
nodes and then over boundary nodes. The associated hat functions are denoted <f>f\ that is <f\ € V) 
and ( t>f{y \) = 1 if 1 = 1, otherwise 4>i(y l f) = 0. Set ( fif := <t>f/\\(l>f Hz, 1 (f2j• Thus, the <f>f are normalised 
in the L°° norm whilst the iff are normalised in the L 1 norm. The mesh size, i.e. the largest 
diameter of an element, is denoted A Xi. It is assumed that Ax t — y 0 as i — > oo. 

Let hi be the (uniform) time step size used in conjunction with Vi, with T/hi € IN, and let sf 
be the A:th time step at the refinement level i. It is assumed that hi —>• 0 as i —>• oo. The set of 
time steps is S) := jsf : k = 0,..., T/hff. Let the f'th entry of diw(sf, •) be 


(diw(sf,-))e 


w(s 


fe+i 




For each a and i, we introduce operators Ef and If to break L a into an explicit and implicit part: 


Ef : H 2 (Q) -»• L 2 (H), w ^ -af Aw + bf ■ Vw + Vfw, 
If : H 2 (Q) L 2 (H), w ^ -Ef Aw + bf ■ Vw + df w, 


with continuous 

A -> C(H) x C(Q, H d ) x (7(0), a ^ (of, bf, cf), 

A C(Q) x C(Q, R d ) x C(Ti), a ^ (Bf, bf, Bf). 

It is required that cf and cf are non-negative and that there is C € R such that 

||cf||L- + ||cf||L- <c, V* € IN, VaG A (5) 


Also, find for each i a non-negative ff which approximates f a : ff ~ /“. The conceptual statements 
L a « Ef + If and f a ps ff are made precise as follows: 

Assumption 1 For all sequences of nodes (yf)i£ w, where in general l = t(i) depends on i: 

“ (df(yf) + af(yf)) || L ao( supp ^) + \\b a - (bf + If ) || £0 o(n j]R d) 

+l|c° - (cf + Bf) || LOO(n) +1 |r - /“iLco (n) ) = °- 

Let (-,-) denote the standard inner product for both of the spaces L 2 ( H) and L 2 (H, H d ), the 
two cases being distinguished by the arguments of the inner product. Consider the following 
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discretisation of Ef and If by operators Ef and If that map H l ( fl) to : for w £ H l { fl), 
£ £ {1 ,... ,Ni = dimV) 0 }, 

(E fw) e := a?(yj)(Vw,V$) + (bf • \7w + *?«;,#>,. (6a) 

{\»t := n?(yf)(Vw,Vtf) + (If ■ Vw + d>,$), (6b) 

(F?)< :=(/“,#>. (6c) 

Throughout this work, we identify Ef and If, when restricted to V), with their matrix represen¬ 
tations with respect to the nodal basis Under this basis, the nodal evaluation operator 

w i->- w{yl) corresponds to the identity matrix Id. 

We now define the numerical scheme for (2). Obtain the numerical solution r>j(T, •) £ Vf by 
nodal interpolation of vt- Then, for each k £ {0,..., T/hi~ 1}, the numerical solution Vi(s^,-) £ Vf 
is defined inductively by 

-dMst •) + sup(Ef Wi (^ +1 , •) + I fvitf, •) - Ff) = 0. (7) 

aeA 

4 Review of monotonicity and uniform convergence 

The proof of gradient convergence in weighted spaces, given in Section 7, is based on the non¬ 
negativity of numerical solutions and uniform convergence to the viscosity solution. 

Assumption 2 For each a £ A, assume that Ef , restricted to If;, has non-positive off-diagonal 
entries. Let hj be small enough so that /qEf — Id is monotone for every a, i.e. so that all entries of 
all h, Ef — Id are non-positive. For each a , suppose that for all v £ V) such that v has a non-positive 
local minimum at the internal node y-. we have (Ifu)^ < 0. 

It was shown in [10, Theorem 3.1] that Assumption 2 implies the existence of a unique numerical 
solution Vi of (7) and that v t is non-negative. 

Let t, = + (1 — d)s-^ +1 £ [sf, sf +1 ] lie between two time steps, D £ [0,1]. Then we interpret 

Vi(t, •) as the linear interpolant between Vi(s f, •) and Vi(s* +1 , •): 

«i(t,0=^i(si,0 + (l-^K(s? +1 ,-). (8) 

Assumption 3 The Hamilton-Jacob i-Bellman problem (2) has a unique viscosity solution v and 

lim ||uj - u|| L oo(n ) = 0. (9) 

In [10] it was demonstrated that Assumption 3, that is uniform convergence, holds if the following 
conditions are satisfied: 

1. Orthogonal projection: Suppose there exist linear mappings Pi : C([0, T], 77 1 (H)) —>• [0,T]xIf 
which satisfy for all (j>f £ Vf 

= Vt G [0,T], (10) 

and there is a constant C > 0 such that for every w £ C' 00 (IL d ) and i £ IN, 

\\PiM\ w ^°°(n) ^ c ||^|lvyi.~(n) and lim \\PiW - u;|liyi>°°(n) = 0. (11) 


4 




2. Boundary control: For each a € A, we define vf: Si —>• V)° to be the numerical solution of the 
linear evolution problem associated to the control a with homogeneous Dirichlet conditions: 
vf(T. •) = Vi(T, •) is the interpolant of vt, and for each k € {0,..., T /hi — 1}, 

(/hi? + ld)«f(s?,-) = -{hi E? - \d)vf{s^\-) + h t Ff. (12) 

Suppose that for each {t, x) € [0, T] x <912 

inf sup lim sup vf (sf , yf ) = 0, 

ae 

where the supremum is taken over the set of all sequences of nodes which converge to (t,x). 

3. Comparison: Let v be a lower semi-continuous supersolution with v(T, •) = vt and u|[o,T]x9fi = 
0. Similarly, let v be an upper semi-continuous subsolution with u|[o,T]x<9n = 0 and v{T, •) = 
vt- Then v<v. 

5 Projection into the approximation space 

For shorthand, let W = W 1,d+1+£ (Qt) 0 L 2 ((0, T), iL?(S2)) C C(Qt) with s > 0. We also use the 
discrete spaces 

Wn= {v € CX^T],!? 0 ) : v\, k fc+ii O is affine in time}, 

which means that functions in Wi have the form of (8) between two time-steps. Observe that 
Wi C W for all i £ IN. 

We introduce the cut-off operation 

C\ : W W, w inax{rt; — ||t> — vf \ L°°(n T )> 0} 

and denote the nodal interpolant on [0, T] x 12 by I*. Finally we define Qi = Ii o C\. Thus Qi is a 
mapping of the type W Wj. Observe that Qj/v € lif satisfies homogeneous boundary conditions. 
Furthermore, from CiV < Vi and the monotonicity of the nodal interpolation operator it follows 
that QiV < Vi- The stability of the max operation and /; gives, cf. [6, Corollary 1.110] with £ = 0 
and p = d + 1 + s, 


IIQHIw < \\ v \\w, , > 

\\Qiv(T,-)\\ L °°(n) <||«(r,.)IUoo(n)- 

Lemma 1 Suppose that v € W. The sequence QiV consists of non-negative functions, satisfying 
homogeneous Dirichlet boundary conditions and Q{V < v t . Moreover, the sequence converges 
strongly in W and L°°(Qt) to v. 

Proof. The convergence of Qj,v to v in W remains. We break the proof into two steps by means of 
the triangle inequality: 


Step 1. Observe that 


II QiV - v\\ w < ||QiV - Iiv\\ w + || Ii.v - v\\w- 

v-CiV = min{||n - vf\ L oo ( n T ),u}. 
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Thus \\v — Civ ||L°o(n T ) —>• 0 as i —>• oo. Moreover, 

||V(u — C'iu)|| i d+i+ e (Q T ) = ||Vv|| L d+i+£( r .), 

where T t = {x £ : u(x) € (0, ||u — i’i||L°°(Q T ))} as we have V(u — C^v) = 0 in Hr \ T*. Because 

n, T* = 0 it follows that ||u — Ci.v\\w —> 0 as i —> oo. Hence, owing to (13), || Q\v — Iiv\\w = 
|| Ii(CiV — n)||u/ vanishes as i —>• oo. 

Step 2. Let 8 > 0. Recalling the density of smooth functions in W there is a vs £ C°°(£It) 
such that ||w ,5 — v\\w < 8. Owing to the uniform continuity of the partial derivatives of vs on the 
compact set fV there is a j £ IN such that 

(maxdfctMf, x)) — (mindfcuHf, x)) < 8 

V X^lT 7 V X(E:T 7 

for all i > j , t € t £ [0, T], and k € {1,..., d}. Let /i = (yu-i,..., fid) € !R d denote vectors with 
11 fi || 2 = 1. Then there is a maximiser x' £ r and a minimiser x w € r such that 

(maxd^u^fjx)) - (minc^u^f, x)) = d^vs{t,x’) - d^vsft, x") 

= ’Y^L L k{d k vs{t,x') - d k vs(t,x")) < y \fi k \8. 
k k 

From the mean value theorem it follows that directional derivatives d^IiVs in the direction of an 
edge of r are attained on that edge by d^vg and thus \\d^IiVg — <9/^,5 ||l°°(t) A ||//||i<5 < yfd8. For 
each t a basis fi 1 ,..., fi d of IR d can be formed of unit vectors parallel to edges. Let A r be the 
matrix expressing the coordinate transformation from fi 1 ,... ,fi d to the canonical basis. Because 
of shape-regularity there is a constant C which bounds the row-sum norm: ||H. r || row < C for all 
t £ % and i £ IN. Thus we have \\d k IiVg — d k vs\\L°°{n T ) < CVd8 for k £ {1,..., d}. We conclude 
that for i sufficiently large 

\\hv — v\\w < || Ii(v — U5)||vy + \\IiV$ — vg\\w + ||nj — v\\w ^ 8, 

using once again (13). □ 

We shall also require a super-approximation result on nodal interpolation in weighted Sobolev 
spaces. We cite Theorem 2.1 in [4], with the notation adapted to this paper. 

Lemma 2 Let T be a mesh element with diameter Ax < 1 and I be the nodal interpolation 
operator of T. Given g £ W 2,00 (T) there exists a value K , depending on ||5 , ||rv 2 .° o (T) and the shape 
regularity of T and the dimension d, such that for affine functions w 

\\g 2 w - I{g 2 w)\\ H i (T) < iFAx(||V( 5 u;)|| L 2 (T) + \\w\\ L 2 {T) ). 

6 Coercivity properties of the Hamilton—Jacobi—Bellman operator 
and its discretisation 

Owing to the non-negativity of v, for each a € A, we formally have for the exact solution 

dtv + sup(L"n — /“) = 0 =► d t v + L a v<f & =^ (d t v, v) + (L & v, v) < (/*, v). (14) 
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Furthermore, if there exists an a € A such that a a € W 2,00 (fl) and c“ — ^(V • b a + A ci a ) > 0 we 
have for w G Hg(fl) 

(L & w , w) = (a & Vw, Vui) + ((c“ - ±(V • b & + Ao“))« 7 , «/), (15) 

thus giving in combination with (14) control on v in a Sobolev space weighted by a a . We intend to 
build the gradient convergence proof upon a bound similar to (14), with the differential operator 
replaced by its discretisation and u by — QiV. 

Fix an arbitrary a € A. It is useful to view Ef and 1“ as bilinear forms on i4 1 (fl) x V). Functions 
u € Vi have the nodal representation 

u(y) = ^2u(j4)</>i(y). 

i 

To test with functions other than <\>\ we introduce the following bilinear form as a partially discrete 
operation: for w € H l (El) and u G V) 

«E fw,u)) :=^2u(yf)(af{yf)('Vw,V(l)i) + (bf • Viu + ^ic,^)). 

e 

We use the corresponding interpretation for ((I fw,u}} and also 

((w,u)) = ((Id w,u)) = Y^w{yl)u(y{)\\4\\^{n), 

l 

e 

Let Hi ( El ) be the closure of C'q°(Q) in the norm 


||u|| 2 := / n 2 7d.x+ / jVnpydx, 

where 7 : fl —>• R is a non-negative L°°(fl) function. We write i^(fl) as abbreviation of H^ a (Ei) 
and Hi (fl) for 77*.(fl) where 7 , is a weight depending on j G If. 

We now formulate a discrete analogue of (14) with (15): Consider that there exists an a € A 
and weights 7 \ and a C' > 0 such that for all i € IN 


w 


£*((0,T),Hl(fi)) 


X 1 


~ Y (((( h i E i - ld )*"( s i + V) + {hi\f + Id )w(sf,-),w(s?, 


k=0 


+ ±((w(T,.),w(T,-)))+C'UVw(T,.)l\ 2 w 


(16) 


(*) 


Y (^« E M^ +1 >•) + 'X4'>■)>^(»i.•)» + H w ( s i +1 ’•) - «>(«*>•) - w(*i> 


k=0 


+ ±((w(0,-),w(0,-)»+C'llVw(T,.)U 2 w 


for all w € IF,; with w > 0 and i € IN, where (*) is a simple reformulation in terms of a telescope 
sum. 
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Example 1 Suppose that there is ana E A such that yfaA G fE 2,00 (n) and c a — |(V-6 Q +Aa“) > 0. 
Choosing a fully implicit scheme with If = L a + 2 1\ Ax* (||-\/a“||vy 1 >«>(r 2 ) + 1) A and E a = 0, the 
highest order term in ((\fw,w)} is at time s\\ 


E i/f)«?(«i, yf)(Vw(st •), v^) 

l 




with re G If 7 ,; as well as weight and numerical diffusion coefficient 

7* : = a? = a" + 2/C Az* (||\/a“||wi,<x»(n) + 1). 
According to Lemma 2, for i sufficiently large, 


(17) 


| (Vw,VIi(a a w)) - (Vw,Va a w)\ < ||Vw;|| L 2 (n) • \\h{a a w) - a a w\\ H i {Q) 

< K A Xi ||Vw|| i 2 (n) (||V(v / a“w)||L2 (n) + |M| L 2 (n) ) (18) 

< K Axi (||v / o“||w r i.<» ( n ) + !) IM|^i (n) . 

Therefore, adding and subtracting (Vio, Va a w) from (17), the triangle inequality and (18) give 
«l fw,w)) = ((Vw,VIi(a Q w) - (Vw,Va Q w)) + (Vro,Va“w) 

+ 2 K A Xi (||\/a"|| l yi,oo(n) + 1) (Vic, Vw) + (b a ■ Vw + c a w, w) 

> (L a w,w) + K A Xi (||\/a“||wn,°o(n) + 1) (Vie, Vic) > ^ \w\ 2 H i^ , 

implying (16) as the reformulation (*) shows. We used here that af — a a is constant and therefore 
unaffected by nodal interpolation. □ 

Due to the definition of the numerical method and the non-negativity of the Vi, if (16) holds 
then 


--1 


v i\vt(Hl) ~ E ~ ld ) u <( s i +1 >0 + {Wi + \d)vi(s!?,-),Vi(s: 


k =0 


+ Uv i (T,-),v i (T,-))) + C'\\v i (T ,-M' 2 


! /II W 


< 


< 


E ((hi??, Vi(st •)» + \{{vi{T, -),Vi(T, •)» + C'\\vi{T, • 


12 
I W 


k =0 


{T ll/f lliffn) + 1) ll v illz,°°([o,T]xn) + C'\\v(T, •) 


lvr> 


(19) 


using stability of nodal interpolation to bound \\vi(T, -)||w- So (16) guarantees stability in the 
L 2 (H') -norm. 


7 Weighted gradient convergence 

Let us assume for a moment that the approximations QiV G Wi satisfy 
f-i 

E «(/hEf - ld)Q iV (^ +1 , •) + (/*!? + ld)Q iV (4, •), (vi - Qiv){s.I •)» ->• 0. 

k =0 
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( 20 ) 



We discuss the validity of (20) in the proof of Theorem 1. 

With = Vi(s%, •) - Qiv(sf, ■), 

| Vi 

£-1 
(16) h i 

< E «(^ E ? - ld K fc+1 + {hi\? + ld)? fc ,C fc » + W T,h '^ T/hl )) + C'\\e /hi \\w 

k =0 

i -1 

= e «(^e“ - id) Wi (^ +i ,.)+ ( hi \? +id) Vi ( S ?, .),£*»+ m T/hi ,e /hi )) 

0 

£-1 

E «(ME? - ld)Q iV ( S ? +1 , •) + (/*!? + ld)Q^(4, + C"||£ T /*||^ 


fc =0 

T 


k=0 


f-1 
(*) ' 


< E «w>£*» - E «(^E?-id)Q iV (s? +1 ,-) + M + id)Q iW («?,■),£*» 


fc =0 


k=0 


+ U(Z T/hi ,Z T/hi )) + c'U T/hi 


\w > 


( 21 ) 


using in (*) the numerical scheme and that, due to the assumptions on the Qi, the sign of Vi — QiV 
is known. Since 


£ «A»F”,e‘>) < ||/”|| 12 £ Ai(|M»?,0 - v(»?,0IU= + MT-) - 


k =0 


k =0 


< 


ll/f IIL 2 (Q) (IK — v llz, 2 (f 2 T ) + ll' u — QM\l*(P. t )), 


the first term in (21) vanishes as i —> oo. The second term vanishes due to (20). For the two last 
terms recall Step 2 in the proof of Lemma 1 and that Vi is the interpolant of v at time T. Hence 


\ v i — v \l 2 ([o,t],h}(Q)) 0 

as i —>• oo. 

Theorem 1 Suppose the value function v belongs to the space W and there is an a € A and 
weights 7 , 7 j € L°°(fl), i € IN, such that 

1 . Ilv'hfllw^n) and \\\/W*\\w 2 ’°°(n) are uniformly bounded in i, 

2 . there is a C > 0 such that af < Cy* and af < C'ji for all i € IN, 

3. the coercivity condition (16) is satisfied for large i and 7 < 7 j. 

Then the numerical solutions converge to the viscosity solution v strongly in L 2 ([0 ,T], H^(Q)). 


9 



Proof. It remains to show (20). The terms connected to the time derivative in (20) vanish in the 
limit as 


E ((QAsi + \-)-Qiv(st-),Z k )) = E 

k =0 k =0 

< II^HIlrqt) ll£ fe llL°°(n T ), 
using the uniform convergence in f k . Recall that ((I fQiv(s k , -),£ fc )) is equal to 

E^ - Qiv)(slyl)(a?(yl)(VQiv(sl-),V<l>i) + (If ■ VQ iV (4, •) + % Qiv(s k , •),<£ 


( 22 ) 


The lower-order terms vanish due to the uniform convergence of v t — QjV to 0 and the bound 

k \ i =a /o „J„k 


sup ||6“ • WQ i v(s i , •) + cf Qiv(s i , -)llLi(n) < °o. 


We note that 


E( v * - Qiv)(sf,yf)5?(yf) (VQM4,-)ytf) = (VQ i v(4,-),VIi(af(v i ~ Qiv))(44), 

l 

so that in (20) the implicit part of the second-order term becomes 
f-i 

h i ,-T 

E h i (VQiv(4, ■), VIi(af{vi - Qiv))(4, •)) = / (JiVQiV, JfWIi(af(vi - Q iV ))) d t, (23) 


k=0 


where maps any w : [0, T] —> L 2 (12;lR rf ) onto the step function (Jiw)\^ s k s k += w(s k , •). Observe 
that JfS/QiV converges strongly in L 2 (£It] R d ). Owing to Lemma 2, we have the chain of inequalities 


l|V/i(afOlli»(n;R-) — \\V( a ?0\\L 2 (n;R d ) + ll^( a fC — Ii( a iO)\\L 2 (n-,R d ) 

< HV(Sf0llL2(n;R<*) + K As i (||V( v ^0llL2 ( n) + UW^m) 

< (llv / “fllvK 1 .°°(n) + K Axj)|| Y / afV^|| L 2( Q;]R d) + (llv^Ilv7i.°°(n) + K Axj)||£|| L 2( n ) 

at a time s k G [0, T). In combination with Assumption 3 as well as (13) and (19) this gives an 
L°°(L 2 ) bound over JfS7Ii(af(vi — Qiv)) in (23). The convergence 

lirn [ (w, J t VI t (af(v t - Qiv)))dt = - lim f (V • w, Jil t {af(v t - Qiv)))dt = 0 
Jo l ^°° Jo 

with test functions w in the dense subset Cq (fiy; U d ) gives weak convergence of VIi(af(vi~QiV )) in 
L 2 (PLt'i K/*), see [13, p. 121]. Combining weak and strong convergence [14, Prop. 21.23], it is ensured 
that (23) converges to 0 as i —>• oo. A similar argument shows that /ij ((E“Qju(s^ +1 , •), £ fc )) 
vanishes in the limit. Therefore we proved (20). □ 

The Sobolev regularity of the value functions is for example discussed in [8] and [12] and [15]. 
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8 The method of artificial diffusion 


We illustrate now a way of choosing the coefficients of Ef and I" in order to satisfy the assumptions 
of the above analysis. 

For all a £ A we need to impose the conditions to ensure uniform convergence. For one a € A 
we wish to enforce terms to guarantee convergence in a Sobolev space with a weight associated to 
L & . 

Remark 1 One could define a set B C A of multiple indices a for which one wishes to derive 
Sobolev norm bounds, which are associated to different weights 7 " and 7 “. The analysis generalises 
directly. 

Given a function g : R —>• and an element K of the mesh 7f we denote 

Ha- := QClkllioo 

3 =1 

If g is elementwise constant then \g\x is simply the Euclidean norm of g on K. Let Ax K denote 
the diameter of K. We assume that the meshes % are strictly acute, cf. [3], in the sense that there 
exists d £ ( 0 , 7 r/ 2 ) such that for all i £ IN: 

• V<ft\ K < - sin(tf) |V $| K |V^U' W ,1 < dim Vi, £ + l, VK £ 7f (24) 

We choose a splitting of the form a° = d Q + a", b a = b a + b a , c a = c a + c", which does not depend 
on i € IN. It is generally necessary to add artificial diffusion to the second-order coefficients, that 
means that we generally need to determine ^-dependent coefficients af > a a and af > a a . It can 
also be necessary to construct cf > c a . For the other lower-order terms we can use the above 
splitting directly: For all i £ IN set 

bf = b a ,bf = h cf = c° and also /f = f a , 

where all terms are in Cffl), a a and a a are non-negative and all c a and c a are non-negative and 
satisfy inequality (5). 

For instance, one could discretise symmetric terms implicitly, i.e. a a = a a and c“ = c a , and 
screw-symmetric terms explicitly, i.e. b a = b a . With this approach the coercivity properties of L a 
are well incorporated. Alternatively, if there is coercivity which is uniform in a with respect to a 
useful weight 7 , then it is interesting to select If = b independently of a because in this case only 
a linear system needs to be solved at each time step while an 0 (A 37 ) time-step may be preserved. 
We also refer to [9] for further illustrations of operator splittings. 

To obtain uniform convergence we select the non-negative parameters uf and vf such that for 
all K that have y\ as vertex: 

(\bf\ K + Ax A -||cf || L oo (A -)) < v? sin(i?) |V$| A vol(R), (25a) 

(\bf\ K + AxkW 11 lo° (a-)) < vf sin(i?) \V%\k vol(iF). (25b) 

For all a £ A \ {d} we choose af and af both in C(R) such that at the nodes 

af{yl)>a a (yf) + vf 1 
af{yl)> a a {yf) + uf J 
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Moreover, we set cf = c a . 

We now turn to convergence in the weighted Sobolev norm associated with the control a. We 
require that VdA, € W 2 ’°°(Q) as well as 

- i(V- b & + Ah") > 0, 

c“ — |(V • 6 " + Aa“) > 0, (27) 

7 := h" - o“ > 0 

with a" < 7 and a“ < 7 . We choose af and af both in C(S1) such that at the nodes 

af{yl)>a & (y{) + vf 1 
af(yf)>a & {yf) + m j 

such that for large i 

Hi > max jp", 2 K A Xi (||y'hf+ ll\/of llw 1 - 00 ^) + 2 ) + hi\\Va a + 6 a |||oo|. (29) 

Notice the recursive nature of the definition of af which appears also on the right-hand side of 
(29). For Ax’j small enough it is clear that Hi can be chosen so that is satisfies (29). Finally, we set 

cf = c" + A A Xi (j| yj$f || w 1 .°°(f 2 ) + + hi He" ||f,oo. 

Sections 8.1 and 8.2 below show that, under suitable time-step restrictions, (26) gives convergence. 

8.1 Verification of uniform convergence 

Suppose that w £ V) has a non-positive local minimum at an interior node yf. It was shown in [10] 
that then 


(E fw)e < 0, (I fw)e < 0, (30) 

and, if uf and uf are chosen optimally, then for K € % 

"?= 0{sup K {\bf\KAxK+ \\cf\\ L ™(K)Ax 2 K }), | 

Df= 0(sup K {\bf\ K Ax K + \\cf\\ L oo^ K )Ax 2 K }). j 

Note that (31) is consistent with Assumption 1. We turn to time step restrictions for semi-implicit 
and explicit methods which give Assumption 2. The non-positivity of the diagonal terms of /ijEf — Id 
expands to 

1 > hi (af (yf)(Vcj>l V#> + (bf • V<f>\ + cf <t>\, #>) 

= hi(o(af Ax]?) + 0(\bf\ K Axfl) + 0(cf)). 

Therefore the time step restriction imposed by L a is hi < inf^(A x 2 K /af(yf)), yf € AT, if there is 
a non-zero a a and i is large. It is hi < inf^ (Axx/|&f (v^Ik) if all af = 0, i € IN, and there are 
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non-zero bf , and is 0(1) if all af and bf vanish. There is no restriction if also all cf are zero. If 
the scheme is fully implicit, there are no time-step restrictions. 

Assumption 3 holds if there is an orthogonal projection, boundary control and comparison, see 
page 4. The former is essentially a quasi-uniformity assumption on the mesh, cf. [5], the latter two 
on the boundary value problem. In fact, the comparison principle is one of the building blocks in 
the theory of viscosity solutions. 


8.2 Verification of Sobolev convergence 

The main step of this section is the proof that (16) is satisfied. With iu k+1 = w(s k+1 , •) and 
w k = w(s k , •), 

«E fw k+ \w k )) = ( Vw k+1 ,X7Ii{afw k )) + ( b & • Ww k+1 + c“ w k+1 , w k+1 + ( w k - w k+1 )), (32) 

using the interpolation operator as in (17). We also find, as in (18), for large i 

\(Vw k+ \VIi(afw k )) - (Vw k+1 ,Vafw k }\ < K A Xi (|| y/0? || wl>00(n) + l) h fc+ V ( n) ll^ll^(n)- 
Therefore, 

(Vw k+1 ,VIi(afw k )) > (Viv k+1 ,afVw k ) + (Vw k+1 , w k+1 Vdf) + (Vw k+1 , ( w k - w k+1 )Vdf) 

- K A Xi (|| + 1 ) \\ w '' ll-ffRQ) • 

Recalling (27) and using that a a — af is constant, we see that 
(Vw k+1 ,w k+1 Vdf) + (b a ■ Vw k+1 + c a w k+1 ,w k+1 ) = {{c a - i(V • ~b & + A d a ))w k+1 , w k+1 ) > 0. 
Now 


,w — w ;| 

~d||2 ll„„fc+l||2 


1 U W k -W k+1 \\ 2 L 2 
,,k +1 


|((Vaf + b a ) ■ Ww k+1 + c a w k+1 
<^(||Va« + 6“||ico ||Va; fc+1 ||| 2 + ||e*||£oo \\w k+1 \\h) + 2k, .. 

<^(||Va Q + b & |||oc ||Vu; fc+1 ||| 2 + \\c a f Loo \\w k+1 \\\ 2 ) + ^-((w k - w k+ \w k - w k+l } 

where we can use Jensen’s inequality as = 1 is a convex combination: 

; [ ^(w k {yi)-w k+1 {yl)) 2 (t)l{ x )& x = ((w k -w k+1 ,w k -w k+1 )). 

Jn ( 




We summarise 




(E> 



(33) 


„ fc +1 


„fc+1' 
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As in Example 1, for large i, 


*w k ,w k )) > {5, & Aw k ,Aw k ) + ((d“ - i(V • b & + A~a & ))w k ,w k ) 


(34) 


+ K A Xi 2 || Jaf || w i,oo + || Ja%\\ w -i,oo + 3 + ^||Vo a + 6 a ||ioc (Viu\ Vw k ) 


+ [-KA Xi 


' a'iWw 1 ’ + l) + ^^ll^lli 00 ) {w k -, wk )- 


We observe that all terms on the right-hand side of (33) can be bounded by corresponding terms 
in (34) of time k or k + 1 —except if w is evaluated at the final time. We choose C such that 


C'\\w T / hi \\ 2 w >\ 


‘|£a + \K AXi (II v'hf 11^/1,00(0) + l )\ wT/h %\ci) + h i \\j i \7w T / hi \\ 2 L2 

+ ^(||Va“ + 6“|||oo ||Vw T / h ‘||2 2 + Halloo \\vJ T / hi \\ 2 L2 ). 

Then we obtain (16) with 7 * = 7 + 


> 


(hi((E?w k+1 + \?w k ,w k )) + i((u ; fc+1 - w k , w k+1 - + ±((w°,w 0 )) + C'\\w T / h % 

k =0 

T_ 

f2(h i ('y i Vw k ,Vw k )^ + ±((w°,w 0 )) > H 2 L 2 mnH i m , 


k =0 


where the last inequality follows because the composite trapezium rule is exact for functions in W t . 
Finally, observe that 7 < 7 j and that there is a C > 0 such that af < C^i and af < C^i for all 
i € IN, as required in the statement of Theorem 1. 
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